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This report covers the final work completed on the research project 
"Radiative Transfer Models for the Earth Radiation Budget and Climatic 
Studies." The work was supported by the NASA/Langley Research Center 
(Experiment Analysis Branch of the Atmospheric Environmental Sciences 
Division) through research grant NSG 1522. The grant was monitored by 
John T. Suttles of the Atmospheric Environmental Sciences Division. 
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FEASIBILITY OF QUASI RANDOM- BAND MODEL IN EVALUATING 
ATMOSPHERIC RADIANCE 


By 

S.N. Tiwari^ and Navin Mirakhur^ 

SUMMARY 

The use of the quasi-random band model in evaluating upwelling atmos- 
pheric radiation is investigated. The spectral transmittance and total 
band absorptance are evaluated for selected molecular bands by using the 
line-by-line model, quasi-random band model, exponential sum fit method, 
and empirical correlations, and these are compared with the available 
experimental results. The atmospheric transmittance and upwelling radiance 
have been calculated by using the line-by-line and quasi-random band models 
and are compared with the results of an existing program called LOIVTRAN. 

The results obtained by the exponential sum fit and empirical relations 
are not in good agreement v>ith experimental results and their use cannot 
be justified for atmospheric studies. The line-by-line model is found to 
be the best model for atmospheric applications, but it is not practical 
because of high computational costs. The results of the quasi-random 
band model compare well with the line-by-line and experimental results. 

The use of the quasi-random band model is recommended for evaluation of 
the atmospheric radiation. 


^Eminent Professor, Department of Mechanical Engineering and Mechanics, 
Old Dominion iJniversity, Norfolk, Virginia 23508. 


’Graduate Research Assistant, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 


INTRODUCTION 


Extensive study of the radiative transfer phenomena in the Earth's 
atmosphe.‘ic system has been carried out in the last two decades (refs. 1 
to 4) . This is important for the understanding of the meteorological 
process on all scales and the spatial variation in surface temperature 
in the Earth atmosphere. Techniques for measuring the Earth's surface 
temperature include airborne instruments and satellite-mounted radiometers. 
In order to understand and interpret the instrument performance and 
readings, it is desirable to develop radiation mocels and numerical tech- 
niques that account for the absorption and attenuation of actual atmos- 
pheric radiation. Development of accurate models for radiative transfer in 
the atmosphere is extremely important for the Earth radiation budget 
studies (refs. 4 to 6). These models have to be used for simulation and 
interpretation of Earth radiation budget measurements as well as for 
retrieval of various surface and the atmospheric parameters from satellite- 
measured radiances (ref, . Since the radiation budget of the planet 
has been identified as an important element of the climate system, its 
measurements are being attempted with increasing accuracy (refs. 8, 9), 

As a result, considerable improvement is warranted in the accuracy of the 
theoretical models dealing with atmospheric radiation transfer. 

Many models for radiation absorption by molecular gases are available 
in the literature. The simplest one is the gray gas model (or the emis- 
sivity approximation] and the most sophisticated and accurate one is the 
line-by-line (LBL) model (or the direct integration procedure]. Between 
the emissivity approximation and direct integration method lie several 
narrow and wide band models and band model correlations which vary greatly 
in complexity and accuracy, A comprehensive review on various line and 
band models is available in reference 4. Use of either an LBL model or a 
narrow band model is suggested for most atmospheric applications. The 
narrow band models usually recommended for atmospheric studies are the 
Elsasser (or regular) model, statistical (Meyer-Goody or Goody) model, 
and quasi-random band (QRB) model. The QRB is probably the best band 
model to represent accurately the absorption of a vibration-rotation band 


and is suitable for calculating the atmospheric transmittance and upwelling 
radiance. The fundamental features of the QRB are discussed, in detail, 
in references 10 to 12, and the procedure for calculating the atmospheric 
transmittance and upwelling radiance is gix'en in reference 12. In spectral 
ranges where both line absoirption and scattering are important, a widely 
used approximation for calculating spectrally integrated radiatii^e flux is 
the exponential-sum fitting of transmissions lESFT) method. The basis for 
this method is that the transmission function for a given spectral interval 
is fit by a sum of exponentials. The method is described in reference 13. 

Radiative transfer models used in earlier climatic investigations 
employed radiation charts, generalized absorption coefficients, and emis- 
sivity approximations (refs. 5, Id-lbl. Rodgers (ref. 17) has indicated 
that the use of multi-interval narrow-band radiative transfer schemes in 
climate modeling studies will constitute a significant step forward and 
result in improved accuracy of the model output. Pels and Kaplan (ref. 18) 
have investigated the effects of using different radiative transfer schemes 
on the thermal structure of the atmosphere and its consequences to atmospheric 
dynamics. They employed two different radiative absorption models, the 
emissivity approximation, and Goody's statistical band formulation, and per- 
formed numerical experiments with the .\’C.'\R general circulation model. They 
observed a significant difference in the cooling rates in the two e.xperiments 
which resulted in significantly different mean temperature fields and 
meridional circulations. 

Very high accuracy can be achieved in the radiation computation by 
using the LBL integration procedure in the radiative transfer models 
(ref. 19). However, the procedure is too cumbersome and makes excessive 
demands on computer time. Tiwari and Gupta (ref. 20) have shoim that the 
QRB model can be used for computing atmospheric transmittances with accuracy 
comparable to that of the LBL method and with computer usage more than an 
order of magnitude smaller. Kunde (ref. 11) has also used this model to 
compute outgoing infrared radiances from planetary atmospheres. However, 
before use of the QRB model can be recommended for Earth radiation budget 
and climate modeling studies, further work needs to be done to validate 
the model on a sound basis. This model should be used for absorption bands 
of different species in different spectral ranges. It is quite possible 


that the model is not justified at shorter wave lengths and smaller pres- 
sure path lengths. Furthermore, under realistic atmospheric conditions, 
the model may give good results in certain spectral ranges but is poor in 
other ranges. 

The objective of this study is to validate the quasi -random band model 
under as many different but realistic conditions as possible by comparing 
the results of this model with available experimental and theoretical 
results. For several molecular species, experimental results for spectral 
tr.ansmittance and total band absorptance are given by Burch et al. (.ref. 211 
under different pressure and path length conditions. Tnus, it is highly desir- 
able to compare the results of the QRB model with these e.xperimental results. 
For cases where experimental results are not available, it is important to 
compare the QRB results with LBL results, For certain spectral ranges and 
atmospheric conditions, results of atmospheric transmittances are available 
in the literature wliich have been obtained by using a sophisticated program 
called L01VTRA.N (ref. 22). It is, therefore, desirable to obtain the QRB 
results e.xactly for these conditions for comp.irison with the L01VTRA.N’ results. 
After these model validations, the aim of this study is to use the QRB model 
for evaluating the upwelling atmospheric radiance under different realistic 
surface and atmospheric conditions. 

The basic formulation of the radiative transfer equations and the e.xpres- 
sions for the upwelling radiance and flux are presented in the next section, 
"Basic Formulation." Details of the spectral models used in the study are 
given under "Evaluation of Transmittance and Integrated Absorption." The 
numerical procedure and data source for calculating the transmittance, 
total absorptance, and upwelling radiance are presented ne.xt i, "Computational 
Procedure and Data Source") . Finally, results of the entire study are 
presented and discussed in "Results and Discussion." 
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BASIC FORMULATION 


Introduction 


The fundamental quantities for describing the upwelling radiance are 
discussed in this section. Following the formulation of the equation of 
radiative transfer, the expression for thermal radiation emerging from a 
plane parallel atmosphere is presented. 

Radiative Transfer Equation 

The amount of radiant energy dE^ , in a specified frequency interval 
(v, V + dv), which is transported across an element of area da and in 
directions confined to an element of solid angle duj during a time dt, 
is expressed in terms of intensity I as: 


dE = I cos -3 dv da du dt 

V V 


( 1 ; 


where 6 is the angle which the direction considered makes with the outward 
normal to da. With reference to figure 1, let us consider a small cylindrical 
element of cross-section da and height dz in the atmopshere. The radiation 
emerges into the atmosphere from the surface whose temperature is at T^ and 
the total thickness of the atmosphere is h. The difference in the radiative 
energy within the frequency interval v and v + dv crossing the two faces 
normally in a time dt and confined to an element of solid angle du) is given 
by 


dl 

V 

ds 


ds dv da dtu dt 


( 2 ) 


where ds is the thickness of the atmosphere in the direction of propagation 
of radiation. This difference in energy must arise from the excess of emis- 
sion over absorption in the frequency interval of the solid angle considered. 

The amount of energy absorbed is equal to 


0 ds"! (I^ dv da dm dt) 


(•') 
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Hero is the density of the absorbing medium and is the frciiuency- 

dependent absorption coefficient. 

The amount of enei’gy emitted is given by 

w' da ds dio dt (.41 

where is the omission coefficient. 

luluating the gain and loss in the pencil of i*adiation, we have 


dl .'ds * - v; I + j v' (.S') 

V V V '^v 

Equation 16 ^ can be rewritten as 

dl /ds * U - I le K' K>) 

V V V V 


where is the nonequilibrium source function. 

.Now consider a medium which is in local thormodiniamic equilibrium. The 
absorption and emission of radiation is given by the Kirchoff's law 


where is the flanck function and is given by 

. 

c‘'[a.\p ihv/cT') - 1] 

Under such conditions, the transfer equation becomes 


171 




dl <ds 

V 


IB, 


- I 

V 


un 


The total flux of radiation arriving at a point in the atmosphere (.from 
above or below point") can be computed by integrating equation (.Ul as shoim 
in reference J3. 


Total Upwelling Radiance 

The radiation emergent from the atmosphere is given by the expression 
(refs. 4, 24): 


E(co) » EqCo)) + Eqj^Coj) + E^(w) + E^^Cu) 


( 10 ) 


where 


E(,(ai) 

ErC«) 


^R*t“) 


thermal radiation emitted by underlying surface and 
atmosphere 

reflected atmospheric radiation from the surface 

incident solar radiation reflected by the surface 

the radiation scattered by single or multiple scat- 
tering processes in the atmosphere without having 
been reflected from the surface 

the scattered energy which has undergone a reflection 
from the surface 


The various components of the upwelling radiation are pictorial ly shown in 
figure 2. 


In the spectral region of infrared measurements, the effect of scattering 
and solar-reflected radiation is usually omitted. Hence, the expression for 
thermal radiation emerging -Trom a plane parallel atmosphere can be written as 


E(to) = Eg (to) + Eqj^(w) = e(u) B(to,Tg) t((o,0) 

+ B[(o,T(z)] [dT(to,z)/dz] dz + p (to) F"(to,T) t(oj,0) (11) 

0 

where e(to) is the surface emittance, B(to,T) is the Planck function, T^ 
is the surface temperature, T(z) is the temperature at altitude z, p (to) 
is the diffuse surface reflectance, F’^u),!) is the downward atmospheric 
radiation, and T(f!,0) and x(to,z) are the transmittances from the top of 
the atmosphere to the surface and the altitude z, respectively. The first 
term on the right-hand side of equation (11) represents the radiation from 
the surface; the second term is the radiation from the atmosphere, and the 
third term represents the reflected component of the downward radiation. 


10 










lAiuativin (.11^ '^an be obtained by intevjratinji equation in a manner 
described in reference The contribution of the reflected atmospheric 
radiation from the surface is usually neglected for surfaces with relatively 
high values of surface emittance and for the spectral regions vvhere the 
downward atmospheric emission is small. 

The contribution of the sunlight reflected from the surface is important 
at shorter wave lengths and is given by the component as 

• 7 [1 - cos e HgUl ll21 

where e is the Sun's zenith angle, [1 - ciuj'i] is the ground reflectance of 
the surface, is the Sun irradiance on the top of the atmosphere, 

X « 1 + f(,0l where f(.61 ■ sec 0 for 0 ^ o0“ and Ch t' for ■» oO'’ 

with Ch 0 denoting Chapman's function, and ti.o 1 » ri.o.01 is the trans- 
mission vertically through the atmosphere. 

For radiation budget and cooling rate calculations, however, the required 
quantity is the flu.x density. Upward flux density can he obtained precisely 
by integrating the upwelling radiance over the zenith angle and the 
azimuth such that 

-1 it/.: 

Fi,k',h) » y* S F.iiu.h') sin cos dt' (.ISI 

0 

Integration of equation (.151 h\ using detailed angular distribution of radiance 
is a tedious problem. However, it is simplified considerably for a plane- 
parallel atmosphere and assuming that the source function in equation (111 
is isotr’opic. It is possible with the above assumption to adapt the two- 
stroam appro.ximation whereby the equations of transfer are reduced to only 
two . 

For the purpose of analysis (i.e., radiation model ingl and measurement 
of outgoing flux, it has been suggested to divide the entire long-wave spectral 
range into the following subregions (.ref. 251 : 

(.a) 0.” to 4 n 

ib> 4 to S u 


1 


(.cl S to 17 u 
(dl 9 to 10 u 
(.el 12 to 18 u 
(.£) 18 to 50 u 

Specific reasons for suggesting this spectral subdivision are given in 
reference 25. Different molecular bands and atmospheric effects are 
encountered in different spectral regions. Solar radiation absorption by 
11^0 and COo occurs in the 0.7- to 4-u region. Absorption by the b.3-y 
H 2 O baud is important in the 4- to 8-u region. The effects of IhO con- 
tinuum bands, .N^O, chlorofluoromethanes and aerosols are important in the 
S- to 12-u region. The contribution of 9.b y 0^ is important in the 9- 
to 10-u region. The combined effect of CO 2 and HoO (but mostly COol is 
important in the 12- to IS-u region, ana the effects of the II 2 O pure 
rotation bands are important in the 18- to 50-u region. 

The greatest problem in computing heat fluxes is the integration of 
equation (111 or (131 over the frequency range of interest. The absorption 
coefficient (and, hence, the transmitt ancel is a highly variable function 
of the frequency, and for accurate work it should be evaluated at small 
frequency intervals. Furthermore, within a band which usually consists of 
thousands of rotational lines, the absorption coefficient at any frequency 
is made up of contributions from many lines. In principle, therefore, 
it is possible to calculate the absorption coefficient with very high 
accuracy by summing the contributions of all intei’vening lines, in practice, 
however, it is a very tedious and time-consuming process. For a wide 
frequency range with several bands, each with 2 large number of lines, 
large amounts of computer resources ai*e required. As such, use of simplified 
but accurate models for spectral absorption is highly desirable. Considerable 
efforts have been spent in the past in devising simplified models to overcome 
the problem of numerical integration ovei t.he complicated line structure of 
the atmospheric spectrum. A complete rsviev^f on different absorption models 
is available in reference 4. The absorption models used in the present 
study ire discussed in the ne.xt section. 
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EVALUATION’ OF TR.\NSMITTANCE .AND INTEGRArhU ABSORPT.ANCE 


Introduction 


.Accurate calculation of the atmospheric transmittance is important for 
the correct evaluation of the atmospheric radiative heat flux. In cases 
where high spectral resolutions are desired, the actual distribution of 
line positions and intensities become important, and the contributions 
fi'om the wings of distant lines need to be included. In such cases, a 
line-by-line model should be employed (ref.>?. 4, li), 201. If the integrated 
signals are measured over relatively wide spectral intervals, an appropriate 
band model could be employed. .A brief description of the spectral models 
used in this study is presented in this section. 

In order to describe the infrared absorption characteristics of a 

radiating molecule, it is necessary to consider the variation of the spectral 

absorption coefficient for a single line. In general, for a single line j 

centered at the wave number a)., this is expressed as 

J \ 

] 

{U) 

1 
I 

where S. is the intensity of the jth spectral line and is given by | 


u)1 


S. f . (u).v .1 


/ 'TO 

I*. . d i t-1 — oJ . 1 

J 


U51 


The quantity f ^ luj - ,o.1 is the line shape factor for the jth spectral 

line. It is a function of the wave number and the half-widtli > . and is 

normalized on u - uj. such that 

J 


r 


f.(u) - uj.l dlul - UJ.I 
J J .1 


lltll 


Several approximate line profiles have been described in the literature. 
The most commonly used profiles are rectangular, triangular Lorentz , Doppler, 
or Voigt icombined Lorentz and DopplerT profiles. The study of line shapes 
and line bx’oadening is an active research field. For various reviews on the 
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p 


subieot, one should refer to references do to 51. Lorento, Doppler, and 
Voigt profiles are of special interest in the atmospheric studies, and 
these are discussed in some detail in reterence 4, The line profile 
usually employed for studies of infrared radiative transfer in the Harth's 
atmosphere is the Lorenta pressure-broadened line shape; the shape factor 
for which is such that the absorption coefficient is -xpressed as 


S. % . 

1 J 


[lu - ••• v-1 


\y'here represents the Lorenta line half-width. From simple kinetic 

theory it may be shown that ■Vj varies with pressure and temperature 
according to the relation 

'j =• \io 

where is the line half-width corresponding to a reference temperature, 

and a pressure F^. The values of m and n depend, in general, on 
the collision parameters and on the nature of the molecules. The values 
of m « 1 and n * 0.5 are employed for most atmospheric studies. 

Hie radiative transmittance at any wave number j is given by the 
relation 


exp(-J'^‘c^ u) 


where c is the absorption coefficient at j per cm-vatm and u represents 
a) 

the pressure path length of the absorber in cm-atm. An appropriate spectral 
model for the absorption coefficient is employed to calculate the transmit- 
tance from equation (191. The total absorption of a single line in an 
infinite spectral interval is given by 

A. = (1 - t 1 d(.jj - uj.) l.dOl 


For a homogenous path, this can be expressed as 
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1211 


A. » P [1 - exp(“ k; . u)] Ju - J.1 
J .j, -*J J 

It should be noted that the integration in equations (.20) and (21 1 extends 
only to the spectral range of a single line. 

The spectral absorptance of a narrow band (consisting of a sufficiently 
large number of spectral lines) may. be expressed by 

A « 1 - T » 1 - exp(l - < dul (22) 

0 

For a homogenous path, the total absorptance of a narrow band is given by 

where limits of integration are considered over the narrow band pass, The 
total band absorptance of a wide band may, in turn, be expressed by 

A = y* fl - exp(- k: u)] d(vij - w^) (241 

where the limits of integration are considered over the entire band pass and 
oi^ is the wave number at the center of the wide band. A brief description 
of the LBL and QRB models and a few relations for the total band absorptance 
are presented in the following subsections. 

Line-by-Line (LBL) .Model 

For a homogeneous path, the monochromatic transmittance as given by 
equation (19) becomes 

T(ii)) = exp(- ic(oi1 u] (251 

The value of the absorption coefficient at the wave number due to direct 
and wing contributions is given by 

>c(y)) = tc^(aj) + v"'(u)) (2(51 
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I 
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D W 

where v u ' and ►- vj' represent the direct and wing contributions resnec- 

tively. The lines whose centers fall close to the wave number under 

cons ivlerat ion are treated as directly contributing lines. The direct con- 
tribution from lines of Lorent , shape is given by equation The wing 

contributions to cul are from lines located far from the wave number under 
consideration. The absorption coeffic'ent in the wing region is calculated 
from the relation 








s 


In the present study, the Lorentc LBL model is used in calculating the trans 
mittance values in certain cases for comparison with the QRB model results. 


Quasi-Random Band iQRBl Model 


The QHB model, introduced by Wyatt et al. (ref. 10), is one of the best 
models available for evaluating the transmittance and absorption of a 
vibration-rotation band. The use of this model in transmittance calculation 
results in a considerable reduction in computational time. The procedure 
for calculating transmittance by the QRB model is discussed here. A listing 
of the computer program is given in Appendi.x A, and various s>Tnbols used in 
the program are e.xplained in Appendix B. 

In this method, the entire frequency range A is divided into a number 
of small subintei'vals S of equal spectral width. The lines in each of 
these smaller intervals are assumed to be arranged in a random manner. The 
average transmittance over the entire spectral range is obtained by averaging 
the transmittances of all subintervals. By considering the intensity of 
the strongest line, five intensity groups are created. Thus, lines whose 
intensities are within 10*^ times the intensity of the strongest line are 
taken into consideration. The expression for the t'^ansmission can be 
written as 
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where «:?Gj,x) is the mass absorption coefficient for the ith absorbing gas, 
is the mass density, ami x is the depth of the level measured from the 
top of the atmosphere. Equation (281 can be changed into an alternate 
form as 
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and P is the partial pressure, o is the mass density of fhe absorber, 

and 0 ^ is the absorber density corresponding to reference condition (STP) . 

The average transmittance due to a single spectral line over the sub- 
interval 5 may be expressed as 


T(a) = r r exp[- S. u f (a),u.)]do). (321 

The average transmittance, x^, due to all the lines in that decade is given 
as 


f^(io) = ®2CP[- Sj u f (a),a)^.)]duj* 


(331 


where N is the number of lines within the decade and is the average 

intensity of all the lines within the decade. The average transmittance due 
to all lines in the five intensity decades cf subinterval 5 is e.xpressed 
by 

T Coj) = f 3xp[- u f (u,u).l]da. I (341 


where k represents the kth subinternal of total span A, Further information 
on this model can be obtained from reference 12. 
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Exponential Sum Fit Method 


Numerous approximate methods have been developed for calculating 
spectrally integrated radiative flux. These methods completely avoid the 
cumbersome procedure of taking care of exact line strengths and line shapes. 
One such method ij the exponential sura fitting of transmissions lESFTl . 

This method which has evolved recently is based on the fitting of the 
transmission function averaged over a spectral interval Au by a sura of 
exponentials (.ref. 13). The ESFT has been used for evaluating the total 
band absorptance of water vapor bands. 

The empirical relations for total band absorptance of 0.94, 1.38, 1.8T, 
2.7, 3.2, and 6.3 u water vapor bands have been obtained by Howard et al. 
(ref. 32) . These relation.^ were obtained by curve fitting the e.xperimental 
results. The e.xperiraental results were represented by two separate relations 
as 
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A =/ A d.j = C + D log u + K log (P + p), A > A^ (36) 

In the above relations, u is the absorber concentration in gm • cm”*', P 
is the total pressure in mm Hg, and p is the partial pressure of H^O in 
mm Hg. The quantities c, k, C, K and D are the empirically determined 
constants, and A^ is the absorption above which the strong band relation 
[eq. 1,36}] becomes applicable. It has been found by Liou and .‘?asamori 
(ref. 33) that the above relations are not continuous when A = A^. They 
have modified these relations, and the absorptivity from the laboratory 
mea.«'.urements has been fitted into a single relation a.*' 


where x = uP^^*^, x = 10”*^'^ 
o 

The values of the constants c, k, C, D, K, A^ and x are given 
by Liou and Sasamori i,re£. 53). The total absorption of each band is fitted 
by a series of exponential functions 
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M 

A » duj[l - ^ w expC- k X) f38j 

where and are the exponential fit parameters, independent of 
atmospheric conditions, and M is the number of divisions in a 
particular band. The exponential fit parameters k^ and w^ have been 
obtained by Stephens (.ref. 34) by fitting the total band absorption of each 
band into the relation of equation (38) . The constant w^ and k^ are 
tabulated in reference 34 for the water vapor bands. 


Empirical Relations for Carbon Dioxide Band Absorptance 

The empirical relations for total band absorption of the carbon dioxide 
bands have been obtained by fitting e.xperimental results in the equation by 
Howard et al. (ref. 35): 
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where u is the CO 2 absorber concentration, P is the total pressure, p 
is the partial pressure of CO 2 , and c is a constant. It has been shown 
by many investigators that the total absorption of the CO 2 band varies as 
u^/2. It has been shown that total absorption is proportional to (P + p)^ 
where k varies from 1/2 to 1/4. The above e.xpression is valid only for 
weak bands. The empirical expression for the total absorption of the strong 
band is given by the following equation 

r A oL = C + D log u + K log (P + p) (40) 

>/ 10 

where G, D, and K are constants. For small values of total band 
absorptance, most of the bands follow the "weak band" relation [eq. (39)], 
and for large values of total band absorptance the "strong band" relation 
[eq. (40)] is followed. Most of the bands have "transition values" of 
total absorption above which equation (40) applies and below which equation 
(39) applies. 

In the present study, the empirical results have been c'^mpared with the 
experimental as well as QRB model results only for the 2.7-u COo band. The 
various constants used in the above relations are given in reference 35, 


Computer Code L0KTR.\N* 


The computer code LOKTRAN has been develop'd to calculate atmospheric 
transmittance over a broad spectral interval, the 0.»5- to 28. S—- region 
iref. IS). This method has been based on an empirical graphical prediction 
scheme. The transmittance is calculated by LBL calculations and averaged 
over a spectral interval of 20 cm”’. This transmission data, over a wide 
spectral range and for a wide range of atmospheric paths, has bc-en calculated 
and presented in different charts in reference 36. The at.tospheric transmit- 
tance for various spectral intervals and atmospheric paths can be predicted 
from these charts. The LOWTRAN is a computer code where all the spectral 
curves, transmittance functions, etc. have been digitalized. This program 
is c versatile and transmittance can be computed for si.x different tv'pes 
of atmosphere. 


COMPUTATIONAL PROCEDURE .\ND DATA SOURCE 


The numerical procedure for evaluating the spectral atmospheric transmit- 
tance and the upwelling radiance and the data source used for the calculations 
are described briefly in this section. The first step involved in accurate 
calculation of the transmittance is the evaluation of the spectral abso-^ption 
coefficient. Use of equation (17) is made in evaluating the absorption co- 
efficient values. 

In calculating the atmospheric transmittance, the atmosphere is divided 
into a number of layers of equal thickness (in the present case, 1 km). For 
the present study, the top of the atmosphere was considered to be 10 km, which 
is approximately the top of the troposphere. The pressure path length is 
given by the expression 

where Q. . is the volume mixing ratio of the ith constituent in the jth layer, 

1. J 

dzj is the thickness of the jth layer, and P^ and T^ are the pressure 
and temperature at the center of the jth layer, respectively. The transmit- 
tance at location z in the atmosphere is given by 

t(o),z) = exp[-J '^ij^“^‘^^ij 

° j i 

Following the procedure for evaluating the atmospheric transmittance, up- 
welling radiance is calculated by dividing the nonhomogenous atmosphere into 
a number of homogenous sublayers. If the gas molecules absorb in a specified 
spectral region Au, then the upwelling radiance is given by 
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Equation (11) is used for calculating the upwelling radiance E in the pre- 
sent study. 

The line parameters needed for this study (position, strength, line, 
width, etc.) were obtained from McClatchey et al, (ref, 36, 37), The 
"McClatchey Tape" is available at the NASA/Langley Research Center. The 


atmospheric temperature and pressure profiles were taken from the U.S. 
Standard Atmosphere iref. 3S1. The concentration distributions in 
the atmosphere for H/O, C0-, .V'O and 0» were taken from McClatchey et al. 
iref. cu'>. The CO.* and X,*0 are assumed to be uniformly mixed in the 
atmosphere, Rotational and vibrational partition functions, required to 
account ;or the temperature dependence of the line strengths, were taken 
from McClatchey et al. iref. 37). 


RESULTS AND DISCUSSION 


Introduction 

Spectral transmittance and total band absorptance of selected molecular 
bands il.vS‘'-u H.'O, CO^' , U.'O, 4.5-u N\'0, and ‘’.bh-u CH^I are 

evaluated by using the QRB model. For a few cases, LBL procedure is also 
used in calculating the transmittance. The results are compared with 
available e.xperimental results. Similar comparisons for .some other molecular 
bands .ire .ivailable in reference dO. The total band .vbsorptance calculated 
by using the liSFT method and eim'irical correlations are compared with the 
e.xperimental as well as QRB model I'esults for the U.'O and CO.' bands. 

Atmospheric transmittances are evaluated for selected spectral r.mges 
of the longwave region by using the QRB model, and these are compared with 
the LOOTRAN results. Finally, the results for upwelling radiance, calculated 
by using the QRB model, ai’e presented to emphasize the use of the QRB model 
for actual atmospheric applications. 

Transmittance and Total Band Absorptance 

The spectral traiusmittance and total band absorptance re.sults are 
presented in this section for different molecular hands. These results 
have been obtained by assuming a homogenous medium with censtant pressure 
path lengths. This is because e.xperimental results were available only for 
homogenous conditions. The results obtained under actual atmospheric 
conditions are presented in the next section. It should also be pointed 
out here that LBL and QRB model results have been obtained for e.xactly 
the same conditions iphysical .ind spectral! for which e.xperimental results 
were available. As such, tiMnsmittance results for certain bauds may not 
be .spectrally symmietric. 

". 00-11 CUii band . - The spectral variation of tlie transmittance, as 
calculated by the QRB and/or LBL models, is illustrated in figures and 4 
for two different pressures but for the same pressure path length. I’lie 
e.xperimental results for exactly the same conditions are also shown in 
these fi.guros for comparison. The overall agreement between these results 
is seen to bo i\uitc good. The transmittance values shown in figui'e 4 loi* 
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Figure 3. Comparison of transmittances of 7.o6-u CH 4 band (.P = 3050 mm 
Hg. u = 2.52 cm-atml . 





Figure 4. Comparison of transimttances of 7.66-u CII 4 band (.? = 34.5 nan 
Ilg, u = Z.SZ om-atml . 


26 


a lower pressure are higher than the values sho\m in figure 3. This is 
especially true in the center portion of the band. This is because, for 
a fi.xed path length, absorption usually increases with the pressure until 
the large pressure limit is achieved. At this point, the central portion 
of the band becomes opaque and radiation absorption occurs only in the 
band wings. At lower pressures, however, absorption of radiation occurs 
in the entire band pass. 

The total C°r integrated) band absorptance results are illustrated in 
figure 5 as a function of pressure path length for three different pressures. 
For lower pressures (P = 50 and 200 mm Hg) , the agreement between the QRB 
and experimental results is seen to be very good for all path lengths. For 
the moderate pressure of P = 300 mm Hg, however, the results are seen to 
be significantly different in the range of intermediate path lengths. It 
is possible that in this range the QRB model accounts for the absorption 
in the band wings more so than the experimental measurement. If results 
were available for sufficiently high pressures, only small differences 
between the two resul'"s would be noticed. This is because at sufficiently 
large pressure the central portion of the band becomes saturated. Although 
it is not convenient to show the limiting (linear and logarithmic) results 
in figure 5, it is noted that all the results tend to approach the linear 
limit toward the small path length and the logarithmic limit toward the 
large path length. Some comparative resulLs for band absorptance are given 
in table 1. These results also e.xhibit the general trend seen in figure 5. 

It is noted from the results of the table at P = 3050 mm Hg that for this 
band the limit of large pressure will be reached at sufficiently high 
pressure. 

4.5-u NnQ band . - The variation of spectral transmittance with wave 
number is sho»m in figures 6 and 7 for two different pressures but the 
same pressure path length. The transmittance results have been calculated 
by the QRB as well as LBL model, and are compared with the experimental 
results. The agreement between the results is seen to be good over the 
entire band pass. 

The band absorptance (i.e., the total band absorptance) results are 
illustrated in figures S(a) and 8(b) for different pressures. Except for 
the results at P = 10 mm Hg, a general agreement between the QRB and 
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Table 1. Comparison of total band absorptance for the 7.66-u CH 4 band. 


p 

(ram Hg) 

u 

(cm-atm) 

(cm” ^ ) 

'^ 9 . 

(cm 

(Aj - A^)/Ag 
(•») 

3054 

2.52 

102.0 

119.2 

-16.7 

301 

2.52 

61.7 

70.4 

-14.1 

55.5 

2.52 

35.0 

37.8 

-8.0 

34.5 

47.30 

115.0 

99.2 

1.3.7 

Ag = / A du) 

(Exp); Aq 

= y* A d(i) 

(QRB) . 



29 



2100 2200 2300 

WAVE NUMBER, cm 


Figure 6. Comparison of transmittances of 4.5-u N-^O band (P = 64.5 mm 
He. u = 5.8 cm-atm) . 
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EXPER IMENTAL (BURCH ET AL . . 1%2 ) 



Figure 8. Comparison of total band absorptance of 4.5-u N.O band. 
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oxperiMcntal results is noted for all pressures at all path lengths. It is 
further noted that the results for P « "oO mm Hg and P « 3000 iran Hg are 
almost identical for the entire pressure path lengths. This indicates tliat 
for this band the limit of large pressure is reached for any pressure higher 
than P 760 mm Hg, For this band, the band absorptance results are also 
given in table 2 fcr certain pressure and path length conditions. In most 
cases, the QRB and experimental results are found to be in very good agreement. 

2."3-u H O band . - Tlie spectral transmittance results of the QRB model 
are compared with the e.xperimental results in figure 9 for P * 862 mm Hg and 
u » 0,101 pr cm. The agreement between the results is seen to be excellent 
over the entire band pass. 

The total band absorptance results calculated by using the QRB model 
and ESFT relation, equation (38), are compared with the experimental results 
in figure 10 for P * 20, 250, and 760 mm Hg. The QRB and ESFT results are 
seen to be in good agreement with the experimental results at higher pres- 
sures (P > 250 mm Hg) . For P * 760 mm Hg, the QRB and ESFT results are 
found to be almost the same as the e.xperimental results. A similar trend is 
also seen from the band absorptance results presented in table 3. As such, 
for this band, use of either the QRB model or the ESFT relation would be 
justified at higher pressures (above 250 mm Hg) . For P = 20 mm Hg, however, 
a considerable difference between the e,xperimental and ESFT results is noted 
for all path lengths, but virtually no difference is seen between the QRB 
and experimental results. The use of the ESFT relation, therefore, is not 
recommended for this band at lower pressures. The QRB model, on the other 
hand, could be employed at all pressure and path length conditions. 

1.87-u HoO band . - The spectral transmittance results, as calculated by 
the ORB model, are compared with e,xperimental results in figures 11 and 12 
for 2 different pressures but for the same pressure path lengths of u = 0.101 
pr cm. The agreement between the results is seen to be very good for the 
entire spectral range of the band. It is noted that, in the central portion 
of the band, the transmittance values are lower for the higher pressure 
(fig. 11), but they are relatively higher for the lower pressure (fig. 12). 
This, however, would be e.xpected because at higher pressures the absorption 
is higher. 


Table 2. Comparison of total band absorptance for the 4.3-.* .N^O band 
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igure 9. Comparison of transmittances of 2.7-u H 2 O band (P = 862 mm Hg, 
u = 0. 101 pr cm) . 







Table 3. Ci^mparifon of total band absorptance for the 2.73-|i ll^O band. 
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Figure 11. Comparison of transmittances of l.S7-u H-.0 band tP = 592 mm Hg, 
u = 0. 101 pr cm) . 
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Comparison of transraittances of 1.87-u H 2 O band (P = 129 mm Hg 
u = 0.101 pr cm) . 




The baud absorptance results, as calculated by the QRB model and ESFT 
relation, are compared with the experimental results in figure 13 for 
P = 40, 250, and 760 mm Hg. The agreement between the experimental and 
ESFT results is seen to be poor for all pressures and path lengths. This 
conclusion may also be derived from the band absorptance results presented 
in table 4. As such, use of the ESFT relation is not recommended for this 
band at atmospheric conditions. The use of the QRB model, however, could 
be made (with reasonable accuracy) for this band also. 

2.7-u COo band . - The transmittance results of the QRB model are compared 
with the experimental results in figures 14 and 15 for two different pressures 
and path lengths. In general, the agreement between the results is good. 

It is obvious from the results of these figures that for higher pressures the 
transmittance values are lower because of higher absorption. 

The band absorption results of the QRB model are compared with the 
experimental results in figure 16 for two different pressures (P = 100 and 
760 mm Hg) . The results are seen to be in excellent agreement for all path 
lengths. The band absorptance results have been calculated by the empirical 
correlations available for this band,. and the results are tabulated in 
table 5 along with other results. From a comparison of the results presented 
in this table, it should be obvious that use of these correlations cannot 
be recommended for most atmospheric applications. 

Atmospheric Transmittanc ; 

The atmospheric transmittance is evaluated by tie QRB model as well as 
by the LOOTRAN Program. The results have been obtained for four different 
spectral ranges, and these are illustrated in figures 17 to 21. For this 
comparative study the top of the atmosphere was considered to be at 10 km, 
and all the results were obtained for the standard atmospheric conditions. 

The results for the spectral range from 3300 to 3700 cm ^ are shown 
in figure 17. The contributions of all the important molecular bands in 
this spectral range (2.73-u H 2 O, 2.7-ji CO 2 , 2.87-u N 2 O, and 2.97-y N 2 O) were 
considered in calculation of the transmittance. The QRB results are seen 
to be in general agreement with the LOWTRAN results. 
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Figure 13. Comparison of total band absorptance of l.S7-u band. 
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Table 4. Comparison of total band absorptance of the 1.87-p H^O band. 
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Figure 14. Comparison of transmittances of 2.7-u CO- band (P = 2065 mm 
He, u = 24,4 cm- atm) . 







Table 5. Comparison of total band absorptance for the 2.7-p CO;* band. 
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Figure 17. Comparison of atmospheric transmittance in the spectral range 
from 3300 to 3700 cm 
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The transmittance results for the spectral range from 2500 to 2800 cm"* 
are shown in figures 18 and 19. The transmittance results as calculated 
by the QRB and LBL models and the L0WTRA.M program by considering only 
the 3.17-u water vapor band are shown in figure IS. The results are seen 
to be in good agreement. The results presented in figure 19 are also 
for the same spectral range as for figure 18, but in this case contri- 
butions of other bands (3.57-y O 3 , 3,85-y CH 4 , and 4.5-y .N 2 O) have been 
included in calculating the atmospheric transmittance. From a comparison 
of results presented in these figures it may be concluded that the QRB 
results axe in good agreement with the -LBL and L0KTRA.\' results and that 
the transmittance (in this spectral range) is mainly due to the 3.17-u 
H 2 O band. 

The results for the spectral range from 1800 to 2000 cm”^ are shown 
in figure 20. The contributions of all the important bands in this 
spectral range (6.27-n H 2 O, 5.2-u CO 2 , and 4.5-u N 2 O) were considered in 
calculation of the transmittance. The QRB results are seen to be in good 
agreement with the LOIVTRAN results. 

The results presented in figure 21 are for the spectral range from 
500 to 800 cm"^. The contributions of all the impoi'tant molecular bands 
(20-u H 2 O, 15-u CO 2 , 14.3-u O 3 , and 17-u N 2 O) wore considered in calculation 
of the transmittance. The QRB results are seen to be in general agreement 
with the LOIVTRAN results. 

Upwelling Atmospheric Radiance 

For the standard atmospheric conditions, the results for upwelling 
radiance have been obtained for two different spectral ranges (2500-2800 
cm”^ and 3300-3700 cm"*) by employing the QRB model. The top of the atmos- 
phere was again taken to be at 10 km. The contribution of the reflected 
component of solar radiation was included in calculation of the upwelling 
radiance, but the contribution of atmospheric radiation reflected from 
the surface was neglected. 

In figure 22, the results are pi:, .ented for the spectral range from 
2500 to 2800 cm"*. The molecular species whose contributions are included 
in this range are H 2 O (3.17 y), O 3 (3.57 u) CH 4 (3.85 y) and N 2 O (4.5 -O • 
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Figure 18. Comparison o£ atmospheric transmittance in the spectral range 
from 2500 to 2800 cm"^ considering water vapor only. 


50 






Figure li?. Comparison of atmospheric 
from 2500 to 2800 cm”^. 
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Figure 20. Comparison of atmospheric transmittance in the spectral range 
from 1800 to 2000 cra"^. 
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Atmospheric transmittance 




Figure 21. Comparison of atmospheric transmittance in the spectral range 
from 500 to 800 cm ^ . 
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Upwelling radiance as a function of surface temperature 
(spectral range from 2500 to 2800 cm 





The figure shows the variation of upwelling radiance with the surface 
temperature for different values of the surface emittance. For a fixed 
surface emittance, the upwelling radiance is seen to increase with 
increasing surface temperature. This is because the surface and atmospheric 
emissions are relatively higher at higher surface temperatures. For a fixed 
surface temperature, the upwelling radiance is seen to increase with 
decreasing surface emittance. This is because, for lower surface emittance, 
the reflected component of the solar radiation is larger, and this makes 
the total upwelling radiance relatively higher. 

The results of upwelling radiance as a function of surface emittance 
are illustrated in figure 23 for different surface temperatures. These 
results were obtained for the spectral range from 3300 to 5700 cm"^ with 
contributions of 2.7-u COo, 2.73-u H 2 O, 2.87-w .M 2 O, and 2.97-u N 2 O bands 
included. The results show the same general trend as seen in figure 22. 

For a fixed surface temperature, it is possible for the upwelling radiance 
to increase with increasing surface emittance if the reflected component 
of solar radiation is not included. 
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CONCLUSIONS 


Spectral and atmospheric transmittances have been calculated by dif- 
ferent procedures to validate the quasi-random band model. The spectral 
transmittance and total band absorptance were evaluated for selected 
molecular bands (7.66-y CH 4 , 4.5-u N^O, 2.73-u H 2 O, 1.87-y HoO, and 2.7-u 
CO 2 ) by using the LBL and QRB models and compared with the experimental 
results. The total band absorptance was evaluated also by using the 
ESFT and empirical correlations for HoO and CO- bands, respectively. 
Atmospheric transmittances were calculated by the LBL and QRB models and 
compared v;ith the results of the LOlVTRAN program. Upwelling atmospheric 
radiances were calculated in two different spectral ranges to demonstrate 
the application of the QRB model. 

For homogeneous conditions, comparison of transmittance and total 
band absorptance results obtained by different spectral models indicate 
that, for all bands, the QRB results are in good agreement with the LBL 
model and experimental results for most pressures and path lengths. The 
LBL model provides the best procedure for evaluating the transmittance, 
but is not practical because of high computational costs. The ESFT 
method has been found applicable only in a few oases. The results obtained 
by using the empirical correlations are not in good agreement with other 
results. Hence, use of the ESFT relation and other empirical correlations 
cannot be justified for atmospheric studies. 

Comparison of atmospheric transmittance, as calculated by the QRB model 
and LOlVTRAN program, indicates that the QRB results are in good agreement 
with the LOlVTRAN results for all spectral ranges considered. In the 
spectral range from 2,500 to 2,800 cm results obtained by considering 
the contribution of only the 3.17-u H-O band indicate that the QRB results 
compare well with LBL and LOlVTRAN results. From these comparisons, it 
may be concluded that the QRB model is computationally fast and accurate 
and, therefore, quite suitable for most atmospheric applications. The 
procedure for calculating the upwelling atmospheric radiance has been 
developed by using the QRB model. For the scandard atmospheric conditions, 
the results for upwelling radiance have been obtained to demonstrate the 
use of the QRB model. The procedure can be used for extensive parametric 


studies (surface as well as atmospheric) of atmospheric radiance. The program 
can be adapted easily for the Earth radiation budget and climate modeling 
studies . 


APPENDIX A 


QUASI-RANDOM COMPUTER PROGRAM TO CALCULATE ATMOSPHERIC 
TRANSMITTANCE AND UPWELLING RADIANCE 


Computer Program QSRBAND 


PROGRAM OSRE-iANCK INPUT, OUTPUT, TAPE7, TAPL8 TAPER, TAPElU 
INTEGER X,G,W 

n I MENS I ON WN ( 5 ) , HS ( 5 ) , OR ( ? ) , HL ( 5(J ) , P'CK ( , .=:•;) , U ) , R ADNC: •: 
D I MENS I ON FR 1 ( 2000 ) , S:A 1 < 20t;H;) ) , AL 1 ( 2000 ) , EL. I k 2000 > , 

/FR2 ( 370 ) , SA2 ( 370 ) , AL2 ( 370 ) , EL2 ( 370 ) , FR3 ( 1 iOO ) , 

/SA3 ( 1 6CiO ) , AL3 ( 1 600 ) , ELS ( 1 600 ) , FR4 ( S!20 ) , SA4 ( 320 ) , AL 4 ( S: 
D I MENS I ON TRH ( 50 ,11), TRC ( 50 ,11), TRN ( 50 ,11), TRHO ( 50 ) , Tl- 
/ TRNO ( 50 ) , OPl- ( 10,4), TRO ( 50 ,11), TROZ ( SO ) , RAlJNCE ( 6 ) , I'EMS ( 
D I MENS I ON QN (10,4), TRA ( 50 , 1 1 ) , TRM ( 1 1 ) , TRT ( SO ) 

COMMlON X 1 ( 26 ) , T1 < 26 ) , X2 ( 21 ) , T2 ( 21 ) , FRB ( 50 ) , FRC ( 50 ) , TRl 
/PREC (10), TEMC ( 1 0 ) , QO ( 1 0 ) , DEL , PNTP , TNTP , TEMR , PR , LC , JD , 
READ 1 2 , FRL , FRU , DEL , PNTP , TNTP , TEMR , RP 1 , RP2 , RP3 , RP4 
PR I NT 1 2 , FRL , FRU , DEL , PNTP , TNTP , TEMR , RP 1 , RP2 , RP3 , RP4 
READ 1 1 , L C , MG , LE 1 , LE2 , LE3 , JD , LE4 
PR I NT 1 1 , LC , NG , LE 1 , LE2 , LE3 , JD , L£4 
READ 1 1 , NG 1 , NG2 , NG3 , N04 
PR I NT 1 1 , NG 1 , NG2 , NG3 , NG4 
READ 1 6 , ( PREC ( L ) , L= 1 , LC ) 

PR I NT 1 6 , ( PREC ( L ) , L= 1 , LC ) 

READ 1 2 , ( TEMC ( L ) , L= 1 , LC ) 

PRINT12, (TEMC(L) ,L=1,LC) 

READ 1 0 , ( ( QN ( L , G ) , L= 1 , LC ) , G= 1 , NG ) 

PR I NT 1 0 , ( ( QN ( L , G ) , L= 1 , LC ) , G= 1 , NG ) 

READ 1 0 , ( ( VPF ( L , G ) , L= I , LC ) , G= 1 , NO ) 

PR I NT 1 0 , ( ( VPF ( L ,'G ) , L= 1 , LC ) , G= 1 , NG ) 

READ 15, (XI (W) , W=l,26) 

PRINT15, (XI (W) , W=l,26) 

READ IS, (T1 (W) , W=l, 26) 

PRINT15. (T1 (W) , W=1 , 26) 

READ 15, (X2(W),W=1,21) 

PRINT15, (X2(W) ,W=1,21 ) 

READ 15, (T2(W) ,W=1,21) 

PRINT15, (T2(W) , W=l,21 ) 

READ 12, EMI, ZEN 
PR I NT 12, EMI , ZEN 
READ 12, (WN(N) ,N=1 ,5) 


, SO ; 


) , SL 4 ( 
\’C0 ( 55 Tj ) 

*■ 6 ) 

i , 11 ) 
/P( 10) 


PRINT12, (WN(N) ,N=1,5) 

READ 1/r- (HS(N) 

PRINT16, (H:s;(N),N=l,5) 

READ 1 2 , ( TEMS ( KK ) , KK> 1 , /:• ) 

PR I NT 1 2 , ( TEMS ( KK ) , KK= 1,6) 

READ (7,20) ( FRl ( X ) , SAl ( X ) , ALl ( X ) , ELK ;< ) , X = 1 , LEI ; 
READ ( S, 20 ) ( FR2 ( X ) , SA2 ( X ) , AL2 ( X ) , EL2 ( X ) , X=1 , LE2 ) 

READ ( •=> , 20 ) ( FR3 ( X ) , S A;S ( X ) , AL3 ( X ) , EL3 (X ) , ,< = 1 . LE 3 ) 

READ (10, 20 ) ( FR4 ( X ) , SA4 ( X ) , AL4 ( X > , EL4 ( X ) , X = 1 , LE4 ) 

10 FORM AT (SFl 0.4) 

11 FORMAT (1615) 

12 FORMAT( IOFS.2) 

15 FORMAT ( lOFS. 4) 

16 F0RMAT(SE10.4) 

20 FORMAT ( 2 ( F 1 0 . 3 , E 1 0 . 3 , F5 . 3 , F 1 0 . 3 , 5X ) ) 

RK= ( FRIJ-FRL > /DEL+0 . 1 
KR=RK 

DELA=0 . 5-s-DEL 

FRB( 1 )=FRL 

DO 100 K=1,KR 

FRB ( K+1 ) =FRB ( K ) +DEL 

1 00 FRO ( K ) =FRB ( K ) +DEL A 
LB=LC+1 

DO 101 L=1,LB 
DO 101 K=1,KR 

101 TRA(K,L)=1. 

IF(N01-1 ) 102, 103, 102 

103 CONTINUE 

DO 104 L=1,LC 
VP(L)=VPF(L, 1 ) 

104 0G(L)=QN(L, 1) 

CALL ASHOK ( FR 1 , SA 1 , EL 1 , LE 1 , AL 1 , RP 1 ) 

DO 105 L=1,LB 
DO 105 K=1,KR 
TRH(K,L)=TRO(K,L) 

1 05 TRA ( K , L ) =TRA ( K , L ) »TRO ( K , L ) 

102 I F ( N02-2 ) 1 06 , 1 07 , 1 06 
107 CONTINUE 

DO 103 L=1,LC 
VP(L)=VPF(L,2) 

109 0G(L)=QN(L,2) 

CALL ASHOK ( FR2 , SA2 , EL2 , LE2 , AL2 , RP2 ) 

DO 109 L=1,LB 
DO 109 K=1,KR 
TRC(K,L)=TRG(K,L) 

1 09 TRA ( K , L ) =TRA ( K , L ) *TRG ( K , L ) 

1 06 I F ( NG3-3 ) 1 1 0 , 1 1 1 , 1 1 0 
111 CONTINUE 


I 


' DO 112 L«1.LC 

VP(D=VFF(L,3) 

. 112 ir!ijiL)=iI'N(L,3) 

CALL ASHOK ( FR3 , S A3 , EL3 , LE3 , AL3 , RP3 ) 

, DO 113 L=1,LB 

I DO 113 K=1,KR 

! TRN(K,L)»TRO(K,L) 

1 1 •: TRA ( K , L ) =TRA ( K , L) *TRO ( K , L ) 

ilO IFvN04-4)114, 115, 114 

115 CONTINUE 

DO 116 L=1,LC 
' VP(L)=VPF(L,4) 

116 00(L)=C!N(L,4) 

CALL ASHOK ( FR4, SA4, EL4 , LE4, AL4, RP4 ) 

' DO 117 L=1,LB 

I DO 117 K=1,KR 

' TRO(K,L)=TRG(K,L) 

, 3 17 tRA ( K , L ) =TRA ( K , L ) *TRG < K , L ) 

114 CONTINUE 

I DO 144 L=1,LB 

I TRM(L)=0. 

> DO 144 K=1,KR 

144 TRM(L)=TRH(L)+TRA(K,L) /KR 
DO 143 K=1,KR 

I TRHO(K)=TRH(K, 1) 

TRCO(K)=TRC(K, 1) 

TRMO(K)=TRN(K, 1 > 

TROZ(K)=TRO(K, 1 ) 

143 TRT(K)=TRA(K, 1 ) 

DO 401 N=l,4 

401 GR(N)=(HS(N+1 )-HS(N) ) / ( WN ( N+1 ) -WN ( N ; ) 

DO 402 K=1,KR 

N=0 

403 N=N+l 

IF ( FRC ( K ) -WN ( N ) ) 402, 404 , 404 

404 I F ( FRC ( K ) -WN ( M+ 1 ) ) 405 , 403 , 403 

405 HL ( K ) =DEL* ( HS ( N > +GR ( N ) ( FRC ( K ) -WM ( N ) U 

402 CONTINUE 

CONSi=lS. 625'J<’1 . E— 07 

' CNST=6 . 625->0 . 3 / 1 . 33 

DO 406 KK= 1 , 6 
DO 406 K=1,KR 
RNUM=DEL*CONS-a-FRC ( K ) **3 
EEX=CNST^FRC(K) 

PCK ( KK , K , 1 ) =RNUM/ ( EXP ( EEX /TENS ( KK ) ^ - 1 . ) 
DO 406 L=1,LC 

^ PCK ( KK , K , L+ 1 ) =RNUM/ ( EXP ( EEX /TENC ( L ) ) - 1 . ) 

406 CONTINUE 

! ZEN= ( ZEN/57 . 2957S ) 

G AM= 1 . + 1 . / COS (ZEN) 


bl 


L'lJ 407 KK»tr6 
f'<ADNC e ( J K ) s»0 . 

DO 407 K=1,KR 

RCOM= < 1 . -EM I) * ( C03 ( ZEN > ) *HL ( \ ) * ( TRT ( K ) ) **OAM 
RADMi: < V K , K ) a£M I +F'i: K ( K > , r , I ) «(-TRT ^ r ) i- 1 . E-0 t-Ri: ijM 
DO 407 L=lsLB 

C OMP=Pi; K ( KK , F : 7 L) * ( TRA ( K , L ) -TRA ( } , L- 1 > ) I . E-07 
RADMC ( K.K I K ) =RADNC ( h K , H ) +COMP 
407 PADNi: E ( KK. ) =R ADNi:E ( KK > +RADNC ( KK , K ) 

PRINT S2 
DO 501 K=1,KR 

!?0 J. PRINT S3, FRB < K+ 1 ) , TRHO ( K > , TRCO ( K 1 , TRNO ( K ) , TROZ ( K ) , TRT Y. ) 
DO f-OS K K=1,6 


PRINT 35,TEMS(KK), 

RADNCE ( KK ) 




FDRMAT</,3X,#FREQ. 

•t- • «T» A I A 1 ^ 

H20TRANS. 

/ \ 

C02TPANS. 

N20TRANS . 

OSTRmN 


/ TOTAL TRAMS. It,/) 

S3 FORMAT ( 2X , F 1 0 . 5 , 4 X , 5F 1 0 . 5 ) 
35 FORMAT ( / , 5X , FS . 2 , 5X , E 1 2 . 5 ) 
STOP 
END 


b 




Subroutine ASHOK for Program QSKBANl' 


• U5F.LU * ii'lS M r v fr » ^ Cu > l_s » • 

IMTSOER 

D I MEN!:. : CN r » < 2000 ) . '= A . 2000 ) i EL 2C-<: 0 . -L . 2'i 
/B10< . 50 ) , L. 1 5 1 50 ) , L I E ■. 50 ) , y'P *. 50 s . I ‘ 5 . 5'! , 
/ BL ( 2i:<00 ) , AO’S 1(5.50), •iUPI’E- 1(5, ’^0 ' 

COiMf'’): < I ( 25 / , T 1 ( 2a. ) , X 2 ( 2 1 ; , T2 ( 2 1 ) , f-'F B < 50 , 
/FREC ( 1 0 ) . TEMC (10), OO a 0 ) , DEL , F NTP , TNTF , T£: 
X=0 

DO 121 K=1,KR 
M*0 

’24. I : X . .'!•£ . -^S ) OOTIi 1 22 
X = X+1 

I F V FR ( X ) . LT . FRB ( K ) ) GOTO 3 D-l 
I F ■; FF ( -X ) . G£. Fr.S ^ - 1 > ) OiO TO 1 2^ 

M=f1+1 
00 TO 1 24- 

327 : F ( K . OS . HR ) GOTO 1 20 

:<=x-i 

123 MP(H:)=M 
1 2 i CONT I NUS 

DO 151 K=1,HR 
DO 1.51 1 = 1,5 
151 :iUri3I( I,K)=0. 

LB=LC-rl 

PI =3. 141592.-..54. 

DELA=0.S-4-DEL 

IJLB=0. 

B3Li=0. 

DO 150 MM=1,LC 
L=LC+1-MM 

C3T= ( SORT ( THMP/TEMC < L) ) ) *PREC ( L ) /PNTP 
PART=VP ( L) ,«• ( ( TEMR/T£H!C ( L) ) :«-.+RP ) *2 . ■•;.•''£+ I *■' 
FAOT= 1 . 43P» ( TEhC ( L ) -TEMR ) / ( TEMC ( L ) <-TEP3R ) 
C0N:5T=0 . 1 «TNTP / PNTP 
PL=C0N3T*PREC ( L ) *00 ( L ) /TEP’O ( L ) 

ALA=0. 

DO 201 X=1,LE 
BL rX ) =C3T*AL ( X ) 

ala=ala-^-bl : X ) 

2<!i I I ( X ) =3A ( X ) *PAPT •*£XP ( EL ( X ) *F Ai.. i ’* 
ALA=ALA/LE 
ULA =PL/(PI*ALA) 

B3U=BSU-s-ULA*ALA 
IJL3 =ULB+OLA 
ALB =S3U/ULS 
RHO =AL3/DELA 
X=0 

DO 203 K=1,KR 

IF 3 H’F ( H ) . EQ . 0 ) iDO^OlO.': 

LI3(H '=X-1 


L IE ( K ) =L I B ( K ) +-MP ( K j -1 
X=LIE(K) 

206 CONTINUE 

DO 207 k=i,k:r 

I F < MP ( K ) . EQ . 0 ) GO TO 2 1 9 

JB=l.IB(K) 

JE=LIE(K) 

BIGI (K)=SI (.JB) 

JC=JB+1 

DO 20S J=JC,JE 

I F ( B 1 0 [ ( K ) . GE . S I < J) ) GO TO 20i3 
BIGT (K)=SI ( J) 

20:5 CONTINUE 

DO 209 1=1,6 
IX=-H-1 

209 BIG( I,K)=BIGI(K)*10.«-»IX 
no 210 1=1,5 

N=0 

SSK I,K)=0 
DO 211 J=.JB, JE 

IF(SI ( J) .GT.BIGI I,K) > GO TO 211 
IF<SI ( J) cLE.BIG(I+l,K) ) GO TO 211 
N=N+1 

SS I ( I , K ) =SS I ( I , K ) +S I ( J ) 

NSI ( I , K ) =N 
211 CONTINUE 

I F ( NS I ; I , K ) . GE . 1 ) GO TO 2 1 0 
NSI ( I ,k:)=i 

210 AVSI ( I , K ) =SSI ( I , K ) /NSI ( I , K ) 

GO TO 207 

219 DO 220 1=1,5 
NSI(I,K:) = 1 

220 AVSni,K)=0 

207 . •ONTINUE 

i-J 310 K=1,KR 
DO 310 1=1,5 

3 1 0 SUMS I ( I , K ) =SUMS I ( I , K ) +AVS I ( I , K > *UL A 
DO 212 K=1,KR 
TRG(K,L)=1. 

DO 213 J=1,KR 
TRD=1. 

jA=IABS(.J-K) 

IF CJA.GT.JD) GO TO 213 
2I=FRC(K)-FRC.( J) 

EPST,=ZI/DELA 
DO 215 1=1,5 
NSJ=NSI ( I, J) 

XI=SLIMSI(I, J) 

RES=0. 

7.F (J.NE.K) GO TO 216 
DO 217 W=l,26 

2=RH0#RH0*X I / ( X 1 ( W ) *X 1( W ) +RHO*RHO ) 
IF < Z . OT . 675 . ) GO TO 200 


64 


Y=EXP(-2) 

GO TO 217 
200 Y=0 

217 RES=RES+Y*T1 (W) 

00 TO 215 

216 DO 210 W=1.21 

Zl=RHO#RHO*Xl'/( (EPSI-X2(W) )*(EPSI-X2<l*J) ) ) 

IF (Z1.0T.675.) 00 TO 300 

Y=EXP(-Z1) 

00 TO 218 
300 Y=0. 

2 1 8 R£S=RESH-Y#T2 ( W ) 

RE8=RE8/6. 

2 1 5 TRD=TRD* ( RES<^ #NS J ) 

2 1 3 TRG < K , L ) =TRO ( K , L) *TRD 
212 CONTINUE 
150 CONTINUE 

DO 320 K=1,KR 
320 TRG(K,LB)=1. 

RETURN 

END 
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AITI-NDIX B 


All, 

AI,2 

AL3, 

AL4 

DEL 


ELI, 

FL 2 

EL3, 

EL4 

EMI 


FRB 


FRC 


FRL 


FRU 


FRl, 

FR2l 

FR3, 

FR4| 

GAM 


GR 


HL 


HS 


JD 


KR 


LC 


LEI, 

LE2 

LE5, 

LE4 

NG 



nxrLANATION OF S\’MBOLS USHD IN COMPUTHR PROGRAN 

Symbols Used in Program QSRBAND 

average line width for the different molecules Il^O, CH 4 , N'^O, 
and O3, respectively. 

width of an interval, cm"‘ 

energies of the lower states for the lines of the .uolecules, 
cm’"’'- 

surface emittance 

wave number at the interval boundaries, cm ’ 
wave number at the intei'val centers, cm ^ 
lower frequency lindt of the range 
upper frequency limit of the range 

wave numbers of the lines of the molecules, cm ^ 

1 + f(. 0 ) where f(9) = sec 6 for ^ 60 and f(9) = Ch(9') for 
> 60; ChCa) is Che Chapman function 

gradients used in calculations of IIL 

computed solar_irradiance at the top of the atmosphere in an 
interval, W cm"^ sr”^ 

tabulated values of the solar irradiance at the top of the 
atmosphere, W cm"^ sr”^ 

number of adjacent intervals on both sides of an interval 
from which the contribution is taken into account 

number of intervals in the band 

number of layers into which atmosphere is divided ( 10 ) 
total number of lines in the spectrum 
total number of gases 


Nlil. .\G2l 
NG3, .\G4j 

PCK 

PNTP 

PROG 

QN 

RADNC 

RADNCli 

RCOM 

RPi, RP2\ 
RP3. RP4) 

SAi, sa:I 
SA5, SA4) 

TPMC 

TEMR 

TOMS 

TNTP 

TRA 

TRCO 

TRHO 

TRNO 

TROZ 

I'RT 

Xl. Ti 
T: 


identifying integers for ihe gases 

Planck's function 

pressure at NTP, bar 

pressure at the different layers, bar 

altitude distribution of the molecules H^O, CH 4 , NNO, and 
Oji, respectively 

radiance in eacli interval 

integrated radiance 

radiance component due to reflected solar radiation 

exponent to account for the temperature dependence of the 
rotational partition function for the molecules 

intensity of the individual lines 

temperature at the center of the layers, K 
reference temperature for the lino parameters, K 
surface temperature, K 
temperature at NTP, K 

combined transmittance of all the interferring gases 

transmittance due to CH 4 only between the top of the atmosphere 
and surface 

transmittance due to H^O only between the top of the atmosphere 
and surface 

transmittance due to N 2 O only between the top of the atmosphere 
and surface 

transmittance due to O 3 only between the top of the atmosphere 
and surface 

combined transmittances of all the interfering gases for each 
interval between the top of the atmosphere and surface 

convolution parameters for integration of direct contribution 

convolution parameters for the integration of wing contribution 


(3 


VPF 

ra 

ALA 

ALB 

AVSI 

BIG 

BIGI 

FACT 

JB 

LIB 

LIE 

MP 

.MSI 

SSI 

SUMS I 


vibrational partition functions for H2O, CH4, N^O, and O3 

wave numbers at which the solar irx'adiance at the top of the 
atmosphere is considered 

Symbols Used in Subroutine ASHOK 

altitude-dependent average width of the lines of a molecule 

altitude-dependent individual width of the lines of a molecule 

avex’age value of intensity for the lines in one decade in an 
interval 

intensity values separating the five decades in each interval 

intensity of the strongest line in an interval 

factor used in the computation of altitude dependence of 
integrated intensity 

number of adjacent intervals from which the wing contribution 
is considered 

serial number of the first line in an interval 

serial number of the last line in an interval 

number of lines in an interval 

number of lines in a decade within an interval 

sum of the intensities of all the lines within one decade of 
an interval 

optical thickness between the top of the atmosphere and the 
altitude under consideration 
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